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Abstract 

We study the 2D scattering of axions from an accelerator like quadrupole magnet using the 
eikonal approximation in order to learn whether or not such a setup could serve as a new possible 
method for detecting axions on terrestrial experiments. The eikonal approximation in 2D is in- 
troduced and explained. We also apply the eikonal approximation to two known cases in order to 
compare it with previous results, obtained using Born's approximation, and discuss its correctness. 
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I. INTRODUCTION 



The question of the scattering of axions from a quadrupole magnetic field has recently 
been proposed as a new possible method for detecting axions on terrestrial experiments such 
as the CAST experiment at CERN. Current axion detection experiments use a constant 
magnetic field, that can be generated by a particle accelerator dipole magnet for example, 
to trace QCD axions. In the case of CAST, these axions are assumed to be emerging 
from the sun with a mean energy which is estimated to be E = 4.2 x 10 3 eV These 
experiments use a standard ID analysis (i.e the inverse Primakoff effect) 2J to evaluate the 
axion-photon reconversion. Therefore, it is an interesting question wether a magnetic field 
which is varying over the scattering region, and thus should be analyzed by a 2D formalism, 
could improve the probability of QCD axion detection in the near future. In this work we 
estimate the cross-section and conversion probability of solar axions that scatter along a long 
accelerator like quadrupole magnet, using a novel 2D scattering method which accounts for 
an axion-photon splitting and presented in the means of an axion-photon duality symmetry. 

In a recent publication, two of us have studied the 2D scattering of axions from a magnetic 
field with cylindrical symmetry j^]. In this work, we expand our study to analyze the 
scattering process when the external magnetic field is generated by a magnetic quadrupole 
using the particle anti-particle 2D scattering formalism, introduced in [dj and reviewed in {3]. 
This new 2D formalism uses a duality symmetry between the axion field and the scattered 
component of the photon to define an axion-photon complex field as 

. = a) 

n 

where 4> is the axion field and A is the z-polarization of the photon. We focus here, as in |3J, 
on the case where an electromagnetic field with propagation along the x and y directions and 
an external magnetic field pointing in the z-direction are present. The magnetic field may 
have an arbitrary space dependence in x and y, but it is assumed to be time independent. 
Hence, for the electric field component in the perpendicular direction to the plane we have 
E z = —d t A. In the case where the scattering process takes place in a medium, the photons 
acquire an effective mass. By matching the photon effective mass to the axion mass, the 
duality symmetry is again discovered [3j. However, for convenience let us neglect the axion 
mass and consider the experiment to be conducted in vacuum so we can write the lagrangian 
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in terms of the new canonical variables \l/ and its charge conjugate \J r * 



£ = d^*d^ - l -p{^*dt<S> - , (2) 

where /3(x, y) = gB(x, y) with B(x, y) being the external magnetic field and is the charge 
conjugation of ty. From this we obtain the equation of motion for \I/ 

d^d^ + ipdtV = . (3) 

We therefore have the magnetic ne.d or 0,2 (the U(i) charge), coop.cd to a charge density. 
Introducing the charge conjugation (C.C) p] (i-e \P — > shows that the free part of the 
action is indeed invariant under C.C. When acting on the free vacuum the A and <fi fields 
give rise to a photon and an axion respectively, but in terms of the particles and antiparticles 
(defined in terms of we see that a photon is an antisymmetric combination of particle 
and antiparticle and an axion a symmetric combination, since 

<j> = -^Or + *) and A = - iT) . (4) 

Hence, the axion is even under charge conjugation, while the photon is odd. These two 
eigenstates of charge conjugation will propagate without mixing as long as no external mag- 
netic field in the perpendicular direction to the eigenstates (i.e axion and photon) spatial 
dependence is applied. The interaction with the external magnetic field is not invariant un- 
der C.C. In fact, under C.C we can see that Si —> —Sj, where Sj = J Cjdxdydt. Therefore, 
these symmetric and antisymmetric combinations, corresponding to axion and photon, will 
not be preserved in the presence of B in the analog QED language, since the "analog exter- 
nal electric potential" breaks the symmetry between particle and antiparticle and therefore 
will not keep in time the symmetric or antisymmetric combinations. Moreover, if the corre- 
sponding external electric potential is taken to be a repulsive potential for particles, it will 
be an attractive potential for antiparticles, so the symmetry breaking is maximal. 



II. THE SCATTERING AMPLITUDE IN A 2D EIKONAL APPROXIMATION 

To apply the results of the previous section to some specific systems with magnetic field, 
we write separately the time and space dependence of the axion-photon field as ^>(f,t) = 
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e luJt ipk{r), which yields 

(-co 2 - V 2 + = . (5) 

In order to develop the space dependent term of the axion-photon field, let us consider a high 
energy, non-relativistic scattering. We assume that the wavelength of the ip field is short, 
i.e kR 3> 1, where R is the length scale of the scattering region and k is the momentum of 
the incoming beam, and that |Vb|/-E <C 1, where |Vo| = gB is the averaged magnitude of 
the potential over the scattering region and E is the energy of the incoming beam. Under 
these assumptions we can address the problem by assuming small scattering angles and 
representing our equation in the integral from of the Lippman-Schwinger equation 

Mr) = e ife> + J d 2 r'G { + \r,r')U(r')Mr') , (6) 
where Gq (r, r') is Green's function given by 

° 2J2nk\r-r'\ (2tt) 2 £ -o J k' 2 - P - le KJ 



where k = u for a massless field, writing the spatial part ot the wave function ipk{f) as 

Mr) = e* r >(r) , (8) 
and substituting into Eq.([6]) yields an equation for </>(r) 6] 



(r) = I + e~ ik -r j d 2 r'G^ ] {r,r')U(r')e ik - r ' <P(r') 

... r r^y§i' rj 



1 + v& I d 2 r> J M ^^ UWm = 1 + Hr) 



(9) 



Choosing k = (k,0) and k! = (k cos(#), k sin(#)), the momentum transfer vector q = k' — k 
is just q = (0, k9) for small angles. Changing integration variables from k' to q in the latter 
equation gives 

Hr) = / d 2 r> / d 2 q - ■ U{r>)${r>) , (10) 

(2n) z J J 2k-q + q 2 -ie 

and since |g| <C 1 we can expand I(r) to a power series in terms of \q\ 2 

1 1 



2k-q + q 2 — ie 2k-q — ie {2k-q — ie) 2 



q 2 + --- (11) 
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As a result of this expansion the transmitted part of the wave function and the scattering 
amplitude can be written as a series as well 

<f>(r) = <f>W + <t>W + ... 
f(k,9) = fU+fW + .... 
Now we turn to calculate I(r) to first order, bearing in mind that we chose the incident 
wave to propagate along the x axis, hence giving 

I (1) (r) = j£p I d\> J d qxdq /^X~ y,)) W)W) = 

= ikJ rfV / d ^ § S^ 5 (y - y'W)m = ( i3) 

= ± / d 2 r'Q(x - x')5(y - 2/)C/(r')0(r') = £ J dx>U(x>, y)(j>(x', y) . 

— oo 

Using the latter result we can evaluate <p(r) and the wave function 



x — 1 dx'U(x' y) 

(r) = 1 + i / dx'U(x', y)<f>(x', y) = , 



i(k-r+± f dx'U(x',y)) 

■0fc(r) ~ e -°° 



(14) 



Since we consider here asymptotic scattering, we need to evaluate Green's function under 
the approximation |r| <C \r'\ and hence 



r 2 — 2r ■ r' + r*' ~ ry 1 — 2— • r' ~ r — r*' • f , (15) 



and then, from Eq. ([6]), the wave function can be written as 

„i(fcr+7r/4) r 

Mr) = e iki ' F + \ d 2 r'e- ik ^U(r')Mr') , (16) 

■v8irkr J 

where fcj = fcx and kf = kf is defined to be the scattered wave. Identifying the scattering 
amplitude from the asymptotic behavior of the wave function 

^ k ( r ) = J-/(0W( fcr+7r / 4 ) , (17) 



we get for the scattering amplitude 



1 f ^ , i f dx"U(x",y') 

f(k, 9) = -== / dVe^/-y y [/(r') e 2t -» . (18) 



III. COMPARISON OF THE EIKONAL APPROXIMATION WITH PREVIOUS 
RESULTS 

In this section we apply the eikonal approximation to two cases, a square well potential 
and a magnetic field with Gaussian ditribution, which were addressed in jsj and compare 
the results obtained in |3( by using Born's approximation with the new method presented 
here. 



A. A Solenoid Magnet 

Let us consider a magnetic field generated by an ideal solenoidal current which is de- 
scribed by a step function realizing a uniform magnetic field pointing in the z direction and 
constrained to a cylindrical region around the origin 

{B z , r < R , 
(19) 
, r>R. 

Thus, the potential associated with the square well is U(x,y) = ugB Q Q(x 2 + y 2 — R 2 ) , where 
Bq is the strength of the magnetic field, g is the coupling constant, and u is the energy of 
the incident wave. Then, we obtain 

1 J nrr, // A igB, 



o/ ; dx"U(x", y') = J ^(x'+ ^R 2 - y' 2 ) . (20) 

Using the current limits on the axion-photon coupling constant (i.e g < 10~ 19 eV _1 [7J), the 
energy of solar axions (4.2 eV) and current limits on terrestrial magnetic fields (< 10 2 T in 
the most extreme cases), the condition <C 1 is obviously satisfied. These values for 

the parameters will be used throughout the rest of this work. 
For this potential the 2D scattering amplitude, Eq. ( ITS]) , is 

/(M ) = ^ J dx'dy'e^-^'e^^^^^ , (21) 

where the integration is preformed over the scattering region. Following the procedure from 
the previous section, we evaluate kf — kt for small scattering angles (i.e q = kf — ki w (0, k6)) 
to get 
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f(k,9) = ^t j dx'e^ BoX ' J dy > e ikey'+i 9 B ^R^ (22) 

In order to calculate the total cross section we shall use the optical theorem in 2D. Since the 
latter equation is continuous in 9, one can take 6 = 0. We can further simplify this integral 
by expanding the exponential to a series in powers of gB$R. For reasonable values for a 
terrestrial length scale of the scattering region and the same values mentioned above for the 
coupling constant g and magnetic field strength B the first order approximations for the 
exponent can indeed be justified. This expansion would be done at the end of the calculation 
in order to have a result which is comparable to the Born approximation calculation, hence 
we might as well do it now. Hence, we obtain 



<f = 2y/fha{f(k, 0)} « \g 2 B* J dx> dy>(x> + v 7 ^ 2 - y' 2 ) = 

= L K ({ + 7W) = <23) 

= 8 (gBpfR 3 
3 2 

To obtain the conversion probability, we calculate the ratio between the number of axions 
arriving at the solenoid and the number of photons produced in the conversion process. The 
number of axions hitting the solenoid is given by multiplying the flux of incoming axions by 
the geometrical cross section of the solenoid, given by gq = 2RL, where L is the solenoid 
length. In order to get the 3D total cross-section (i.e the scattering cross-section) a,s we 
multiply the 2D cross-section a tot by the length of the solenoid L. The number of produced 
photons is found by multiplying the scattering cross section (as = o tot ■ L) times the flux. 
Thus, the conversion probability is given by 

P W en = os/a G = l 9 -^f = \g-BlR- . (24) 

Comparing this result to the result obtained by using the Born approximation in j^] for 
the same setup (Eq. (3.20) there) we get precisely 

peikonal i pBorn i (OK\ 

r well I r well — 1 • \ Z0 ) 
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Hence, there is a complete correspondence between the eikonal approximation and the Born 
approximation in this case. 



B. Gaussian Magnetic Field 

In this setup the potential has the form 

U(x,y) =ugB e~ £ W 1 . (26) 
Integrating the potential along the axis of the incident wave yields 



2k fjx"U{x",y>) =;^^[l+Erf(^)]e-£, (27) 



2k J 

where Erf(x) = ^= e t2 dt. The scattering amplitude is given by 



_ _y' 2 ,2 



(28) 

In order to calculate the total cross section we use, as usual, the optical theorem in 2D and 
by using the same reasoning as in Eq. ( 1221) to consider only the 6 = angle we get 



Gauss 



4 / dy' sin 2 {-gB R^e~^) . (29) 

J — oo 

To obtain an analytic result and simplify the calculation, we can use the fact that 
exp{ — } < 1 for all y and that reasonable values of the parameters g, B and R allow us 
to make a first order approximation. Thus, Eq. ( 1291) can be written as 



a 1 



a Gauss „ <_ g 2 BlR 2 I = « ^gtf . (30) 

J — oo \ oZ 

Hence, using the same method we used in the previous setup, the probability will be 

3 

"^222 / \ 

PGauss = ~^~f^9 B Q R . (31) 
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The comparison to the probability of conversion calculated in |3J, using the Born approxi- 
mation, (Eq. (3.13) there) gives 

ryeikonal I ryBorn -r {QO\ 

r Gauss I r Gauss ~ L : \ 6Z ) 

so that there is again a complete correspondence between the eikonal approximation and 
the Born approximation. Thus, we conclude that the Eikonal approximation is indeed valid 
and will most probably yield a correct results under our assumptions for a high energy, yet 
non-relativistic scattering. 



IV. AXIONS SCATTERING IN A QUADRUPOLE MAGNET 

After verifying the accuracy of the eikonal approximation for two known problems, we 
now turn to calculate the scattering of axions from a quadrupole magnetic field. 

Placing the quadrupole magnet in the yz plane with the x axis (the direction of the 
incoming beam) along the symmetry axis of the quadrupole field, the quadrupole magnetic 
field distribution can be approximated (for a quadrupole magnet with a narrow aperture 
compared to its length) by (9j] 

B(x, y, z) = sV(yz) = s(zy + yz) = —j^-{zy + yz) , (33) 

where s = ^ is the quadrupole gradient, B is the value of the magnetic field at the pole 
tips (i.e at the points (0, ±-|, 0) and (0, 0, ±f )), where we chose a rectangular aperture for 
simplicity. However, since in our formalism we chose to analyze the case where the magnetic 
field is pointing in the z direction, only the z component of the external magnetic field will 
take part in the scattering process and we effectively have an inhomogeneous magnetic field 
of the form B z = ^-y- However, the y component of the magnetic field will, of course, give 
the same contribution to the scattering process as the z component with the sole difference 
that the final photons will be with a y-polarization. Thus, we need to take into account the 
^ particle that has the y-polarization of the vector potential as one of its conjugate fields. 
This will be done at the end of this section. 
Defining, in this case, the potential as: 

U = ug^y , (34) 
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for -R/2 <y< R/2 and -L/2 <x< L/2, we get 



± £ dx"U(x", y') = zg^y'(x' + L/2) . (35) 
Putting this into Eq. (118p we have for the scattering amplitude 



f(6) = y/% ^ j^dy'^y'y'e^y' j^dx'B^^ 



2y/Sf% dy'e^^y' sin(^) = (36 ) 



— . f 'sin(|fei^) R Sm(^(gB L+kRd))\ 
- * ' 2tt ' * ' I kd kRO+gBoL 



Therefore, using the optical theorem, we obtain the total cross section 

/ 87F t f^i-oofi sin(^5 L) 




^- =UTM/ (0) } = 2fl^l-^|^j . (37) 
To get the probability, in this case we just divide by R (geometrical cross section = R 2 ) 

sin(k5 £) N \{gB Lf lg 2 B 2 L 2 , N 

Pauad = 2 1 1 2 » 2 8 ^ ^ = ^_ . 38 

Although, as was discussed before, this result may be different than the one that would be 
obtained by using Born's approximation (since the potential is not piecewise continuous), 
it coincides with a result that is obtained by using an optical analogue as was previously 
shown by J. Redondo [10] . However, Redondo's calculation is missing features of the 2D 
calculation (and is also computed in an unphysical setup), as we now show (and as will be 
discussed in the conclusions of this work). 

The scattering from the z component of the magnetic was merely a matter of choice 
in the in the initial setup of our system. The ip field will scatter, of course, from the 
y component as well since this component of the magnetic is also perpendicular to the 
momentum of the incoming beam. The same process described above can be repeated by 
a 7r/2 rotation of the system in the yz plane. In this case, the tp field would have been 
defined as ip = ((f) + iA) / a/2, where A is the ^/-polarization of the photon this time. This 
will give the same expression for the cross-section and conversion probability and since these 
two processes are distinguishable we can sum incoherently the two probabilities. Hence, in 
order to get the complete probability we have to multiply Eq. fl38|) by a factor of two and 
thus 
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p 



V. CONCLUSIONS 



In this paper we have calculated the axion-photon conversion probability from a 
quadrupole magnetic field. We have used the eikonal approximation to calculate the scat- 
tering from the quadrupole magnet which simplifies the calculations compared to the Born 
approximation and verified this approximation by comparing to known results obtained with 
the Born approximation. The comparison to the step function and Gaussian distributed 
fields shows that the eikonal approximation and the Born's approximation give reasonably 
close answers. 



In a previous letter by J. Redondo 10], a related analysis concerning the evolution of 
an axion-photon field in the presence of a magnetic field with a constant gradient over the 
entire space, was studied. In our case, however, the constant gradient field exists only 
in a finite region of space. This boundary condition for the magnetic field contributes in 
an essential way to the scattering amplitude since, as was explained in even for the 
case of a strictly constant magnetic field living in a finite region of space one gets a non 
trivial scattering amplitude due to these boundary conditions. Thus, we conclude that 
the boundary effects can be as important as the gradient of magnetic field inside the finite 
scattering region. This is a significant difference between our treatment and that of reference 
jlo| . In Redondo's letter, the applied magnetic field does not satisfy Maxwell's equations 
without sources since, in his letter, VxB = d x B y = B\ = J z (where in Redonodo's notations 
By = Bix). Therefore, in his analysis there is an infinite extent of currents along the z axis 
which can simply not represent a physical situation. In conclusion, the work of Redondo 
is not a scattering problem and therefore cannot incorporate all the physical aspects of an 
experimental set-up as opposed to the work presented here. This is most easily seen by 
comparing Eq. ( 139|) with equation 3 in Redondo's letter. This comparison unveils the fact 
that the different approaches lead to different results. 

Moreover, it may sometimes be tempting to believe that the 2D results can be obtained 
by averaging over ID conversion probabilities. As was already shown in [3j, a general 
prescription to find a ID analogue to the 2D calculation may be obtained by using the 
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magnetic flux as the averaging measure 



p av 9 = IZo IZo W\ IZo gO^> y)dx'\ 2 B(x, y)dxdy 
1D IT (X> ST o0 B{x,y)dxdy 
However, since the quadrupole magnetic flux is zero, this method will not work this time. 
Of course, it is possible to find an averaging process that will produce the 2D result, but 
it cannot be done a-priori with certainty. Eq. f )40|) is an example for a legitimate choice of 
measure that cannot produce a result at all. It is clear that there is no way to avoid the 
real calculation and obtain results in higher dimensions from averaging on lower dimension 
estimations. This shows that 2D processes cannot be reduced to a ID calculation. The 
scattering process from a quadrupole field is intrinsically 2D since the scattered photon 
may have two different polarizations: Considering a magnetic field produced by a physical 
localized current naturally makes a significant difference. In particular, it shows that any 
attempt to say that the problem can be deduced from a ID analogue is untenable and, in fact, 
will fail since the magnetic field produced by a physical current will necessarily have at least 
two components (in the source free region). Hence, in the process of photon production 
from the scattering of axions, the photon will have two distinguished polarizations. In 
addition, in the source free regions, a non-uniform field pointing in one direction, of the 
form B = B(x, y)l, will not be able to satisfy the source free Maxwell's equations, Vx5 = 
VB(i, y) x I = and V • B = VB(i, y) 4 = 0, since the solution requires that VB(x, y) = 0. 
Therefore, a one directional magnetic field cannot solve these equations and the problem 
cannot be reduced to ID. 

We can, however, compare the scattering from a quadrupole to the scattering from a 
solenoid, which can represent a dipole accelerator magnet (with a different geometry of 
course). The result obtained here shows that it will be preferable to have a constant magnetic 
field distributed over the scattering region aperture (like, for example, the field of a solenoid) 
rather than having an inhomogenous field (as the quadrupole field). This comes from the 
fact that the magnetic field energy will be higher in the first case (when the maximal 
magnetic field strength B is equal for both cases and both fields are distributed over the 
same scattering region). Since B is an intrinsic property of the superconducting material, 
comparing a quadrupole magnet to a dipole magnet of the same length and aperture and 
made with the same superconductor, the dipole will yield a higher conversion rate. 
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One can also observe that the magnetic field parameter that determines the conversion 
rate is actually the global maximum value of the magnetic field in the scattering region. 
This feature appears as well in other calculations of the conversion probability from inho- 
mogeneous fields, like, for example, the Gaussian distributed field in j^]. 
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